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Quantum  Corrections  for  Inelastic 
Atom-Surface  Scattering 

J.  H.  Jensen  and  W.  Kohn 

Department  of  Physics 
University  of  California 
Santa  Barbara  CA  93106 


ABSTRACT:  We  derive  simple  and  exact  expressions  for  the  leading  quantum 
corrections  to  the  angular  and  energy  distributions  of  an  atom  scattered  inelasti- 
cally  from  a  solid  crystalline  surface.  The  quantum  corrections  are  given  in  terms 
of  classical  trajectories  and  equilibrium  correlation  functions  for  the  solid  surface. 
Our  results  are  useful  when  the  incident  atom  has  an  energy  of  ~  1  electron  volt  and 
the  surface  has  a  temperature  less  than  or  comparable  to  the  Debye  temperature. 
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I.  Introduction 


A  neutral  atom  with  an  energy  <  1  electron  volt  colliding  with  a  solid  surface  target  ex¬ 
changes  energy  with  the  target  through  emission  and  absorption  of  phonons. 1  Rapidly  im¬ 
proving  experimental  techniques  have  made  a  detailed  theoretical  understanding  of  atom- 
surface  scattering  a  matter  of  practical  importance.2  Atom-surface  scattering  is  also  of 
inherent  theoretical  interest  as  a  challenging  many-body  quantum  dynamics  problem.3 

Incident  atoms  (“incidons”)  with  small  initial  energies  (~  0.01  ev)  usually  excite  only  a 
few  phonons,  and  the  atom-surface  interaction  is  highly  quantum  mechanical.  With  larger 
initial  energies,  more  phonons  may  be  excited  and  quantum  corrections  to  a  classical 
description  of  the  scattering  process  tend  to  be  smaller.  For  initial  energies  of  —  1  ev, 
quantum  corrections  are  frequently  small,  but  significant,  and  a  semiclassical  approach  is 
appropriate.  This  last  regime  is  the  subject  of  this  paper. 

In  a  previous  publication,4  (hereafter  called  I)  we  considered  the  problem  of  a  particle 
scattered  inelastically  from  a  target  with  internal  degrees  of  freedom.  We  assumed  the 
target  to  be  initially  in  its  ground  state  or  at  a  low  temperature  and  the  particle-target 
system  to  have  a  well-defined  classical  limit.  We  gave  an  exact  method  for  obtaining  the 
leading  quantum  correction  for  physical  quantities  possessing  an  ^-expansion.  The  present 
paper  applies  this  method  to  obik’n  explicit  expressions  for  O(h)  corrections  to  the  angular 
and  energy  probability  distributions  of  an  incidon  scattered  from  a  crystalline  surface.  The 
expressions  are  simple  and  numerically  tractable. 

While  quantum  corrections  can,  in  principle,  be  obtained  with  standard  semiclassical 
methods,3,5  explicit  expressions  for  the  leading  quantum  corrections  have  not  been  previ¬ 
ously  given.  The  purpose  of  this  paper  is  to  derive  such  expressions  and  to  show  that  they 
are  simple  enough  to  be  applied  to  realistic  scattering  problems.  The  main  body  of  this 
paper  assumes  that  the  target  is  at  zero  temperature.  The  straightforward  generalization 
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to  finite  temperatures  is  presented  in  Appendix  A. 

II.  Statement  of  Problem 

This  paper  develops  a  theory  of  quantum  corrections  for  atom-surface  scattering  in 
regimes  where  a  strictly  classical  theory  (Newton’s  equations)  is  a  good,  zeroth  order 
approximation.  Hence  we  begin  with  a  brief  discussion  of  relevant  features  of  classical 
atom-surface  scattering. 

We  establish  a  Cartesian  reference  frame  where  the  origin  is  taken  at  a  point  in  the 
surface  region  and  where  the  z-axis  is  perpendicular  to  the  surface.  We  introduce  spherical 
coordinates  in  the  conventional  way, 

z  =  r  cos  8,  x  =  r  sin  9  cos  <£,  y  =  r  sin  9  sin  4>. 

The  initial  state  of  the  incidon  is  completely  specified  by  six  parameters,  its  initial 
coordinates,  (xq,j/o,zq),  and  momenta,  (pxOtPt/0>Pzo)  (8ee  Fig.  1).  The  “impact  point” 
A  on  the  zy-plane  is  the  point  where  the  incidon  would  penetrate  this  plane  if  it  were 
undisturbed  by  the  target.  From  Fig.  1,  it  is  evident  that  identical  scattering  will  take 
place  whether  the  incidon  is  initially  at  point  B  or  B' ,  leaving  five  quantities  to  specify 
the  initial  state  of  the  incidon.  We  take  these  to  be  the  spherical  angles  (0q>  4>q)  of  the 
incident  momentum  vector,  the  incident  energy,  Eq,  and  the  impact  parameters,  a  and  b. 

A  parallel  beam  of  incidons  averages  uniformly  over  a  and  6.  For  an  incident  beam 
of  given  energy  and  direction  we  denote  by  1(9, 0;  E)  the  intensity  of  the  corresponding 
scattered  beam.  1(9,  <f>;  E)  (sin  9)  d9d<j>dE  represents  the  number  of  particles  leaving  the 
surface  per  unit  time  through  the  solid  angle  element  dfl  =  (sin  9)d9d<£  and  in  the  energy 
range  dE,  divided  by  the  total  incoming  flux. 

A  particle  with  specific  impact  parameters,  (fli,&i),  will  be  scattered  in  a  unique 
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direction,  (0i,0i),  with  a  unique  energy,  E\,  smaller  than  Eq.  However,  there  are,  in 
general,  other  impact  parameters,  (a*, 6*),  which  lead  to  scattering  with  the  same  (0i,<£i). 
Some  of  these  impact  parameters  are  equivalent  to  (ai,6})  and  hence  lead  to  the  same 
scattered  energy,  E i,  while  others  are  physically  inequivalent  with  E^  #  E\  (see  Fig.  2). 
Hence,  for  a  given  8  and  <£,  1(8,  fa  E)  consists  of  a  sum  of  delta  function  peaks  in  E ,  as 
shown  in  Fig.  3.  We  call  /*(£,  fa  E)  the  intensity  associated  with  the  Jfcth  energy  peak  so 
that 

I(»,*-,E)  =  (1) 

k 

We  also  introduce  the  integrated  angular  distribution  for  the  fcth  peak 

4(M)  =  J  It{»,*-,E)dE,  (2) 

and  the  energy  probability  distrib  ution  for  given  for  given  8  and  <f> 

pm- 

*'  (3) 

=S(E  -  Ek). 

The  angular  distribution,  1^(8,  fa),  corresponding  to  a  uniform  distribution  of  impact  pa¬ 
rameters  is  a  smooth  function  of  8  and  <j>,  disregarding  rainbow  scattering  singularities. 

What  are  the  qualitative  changes  in  the  form  of  1(8,  fa  E)  when  quantum  effects  are  in¬ 
cluded?  We  treat  h  as  formally  small  so  that  hup  <  Eq,  where  up  is  the  Debye  frequency. 
An  incident  plane  wave  then  excites  with  finite  probability,  0,  1,  2,  . . .  phonons,  a  typical 
number  being  n  ~  (Eo  —  E^/hup  >  1.  The  total  angular  distribution,  /  1(8,  fa  E)dE,  is 
now  no  longer  smooth,  but  has  delta  function  diffraction  peaks  due  to  elastic  (0-phonon) 
scattering.  These  peaks  have  weights  that  are  exponentially  small  in  h  and  are  separated 
by  small  angles  A0,A0  ~  h/pQd,  where  d  is  a  lattice  parameter  and  pq  is  the  incidon’s 
initial  momentum.  Since  the  frequency  spectrum  of  the  target  is  continuous,  the  energy 
distribution,  in  any  particular  direction,  is  no  longer  a  sum  of  delta  functions,  but  rather  a 
sum  of  peaks  with  finite  widths  near  the  classical  peaks  (see  Fig.  3).  As  we  shall  see,  the 
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peaks  have  widths  of  O(Aj)  and  are  displaced  from  the  classical  positions  by  an  amount 
of  O(h).  The  overlap  of  different  peaks  is  exponentially  small  in  K  and  immaterial  for  our 
purposes. 


Because  of  the  delta  function  peaks  (both  classical  and  quantum),  a  straightforward 
expansion  of  Ik(6,  <£)  and  Pk{E)  in  powers  of  K  is  not  possible.  However,  one  can  define 
averaged  quantities  which  do  have  such  expansions.6  We  define  the  smoothed  intensity 


where  the  angle  a  sets  the  “resolution.”  We  also  define  a  smoothed  angular  distribution 
by 

7*(M)  =  J  lk(0,<t>-,E)dE,  (5) 

and  a  smoothed  energy  distribution 

fIt(»,*;E)dE' 

The  smoothed  angular  distribution  has  an  h-expansion 

=  f V/mIM).  (7) 

n=0 

where  is  the  classical  limit  of  Ik(0,4>)f  hlk  i  is  the  first  quantum  correction  to 

7*,  2  is  the  second  correction,  etc.  Ikn  tends  to  a  definite  limit,  Ik  n,  as  a  — ►  0.  We 

call  h/jfc  i  =  A/jfc  the  leading  quantum  correction  to  Ik{0,4).  It  neglects  the  exponentially 
small  and  densely  spaced  diffraction  peaks. 


PkW 


In  order  to  characterize  Pk{E ),  we  define  the  energy  moments 

=  f  dE  Pk{E)  ( E)n ,  (8) 

Each  moment  has  a  power  series  expansion  in  h 

oo 

Mn  -Y.  Mn, (9) 

n'=0 
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As  the  resolution  a  approaches  zero,  the  quantities  Mn  ni  approach  finite  limits,  A/n  n<. 
Mn  o  i*  simply  (U*)n.  hM\  \  =  A is  the  leading  quantum  correction  to  the  center  of 
the  peak,  and  Afcf2,l  —  i  =  S^E^  is  the  square  width  of  the  peak  to  0(h).7 

The  purpose  of  this  paper  is  to  give  explicit  expressions  for  A /*,  A E^,  and  S^E^. 

III.  Semiclassical  Formalism 

We  now  review  the  semiclassical  formalism  presented  in  I.  In  the  Heisenberg  represen¬ 
tation,  a  quantum  dynamics  problem  requires  solving  for  the  time  evolution  of  quantum 
operators.  In  I  we  showed  that  quantum  corrections  of  0(h)  can  be  calculated  entirely 
from  classical  trajectories  and  are  due  to  the  equilibrium  zero  point  oscillations  of  the 
target. 

This  result  can  be  established  by  representing  quantum  mechanical  operators  with 
so-called  “symbols,”8  which  are  functions  on  classical  phase  space  that  contain  the  same 
information  as  the  operators  they  represent.  Among  the  many  ways  of  defining  such 
symbols,  a  convenient  choice  is  the  Wigner  form 

o*(z,p;t)  =  J  da  e“»Qp  <  x  +  ||A(t)|x  -  |  >,  (10) 

where  A{t)  is  a  time-dependent  quantum  operator  evolving  according  to  the  Heisenberg 
equation  of  motion,  |x  >  is  a  position  eigenstate,  and  o*(x,p;f)  is  the  symbol  for  A(f).9 
In  Eq.  (10),  x  and  p  are  a  shorthand  for  all  the  position  and  momentum  variables  of  the 
system.  Given  o*(x,p;t)  one  can  reconstruct  the  operator  A(t)  from 

A(t)  =  J  dxdp  a*(x,p;t)  j  .‘MiTO-.HntiHO)-,],  (u) 

where  x(0)  and  p(0)  Me  the  position  and  momentum  operators  evaluated  at  the  initial 
time,  t  =  0,  and  N  is  the  number  of  position  variables. 
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The  expectation  value  of  A  is  simply 

<  ^(0  >=  J  dxdp  /0(x,p)  a*(x,p;t),  (12) 

where  /o(x,p)  is  the  Wigner  function10  for  the  initial  (normalized)  state,  |0O  >,  given  by 

/o(*iP)  =  J  <  *  +  flV’O  ><  iM*  “  |  >  •  (13) 

Clearly,  both  fq  and  a*  contribute  to  the  ft  dependence  of  <  A(t)  >. 

An  important  property  of  a*(x, p; <)  is  that  in  the  small  ft  limit  it  approaches  the 
classical  function  oc(x,p;t),  which  is  the  value  of  the  classical  quantity  corresponding  to 
A(t),  for  initial  conditions  given  by  the  positions  x  and  momenta  p.  For  example,  if 
A(t)  =  x(t),  then  ac(x,p;t)  =  xc(x,p;t)  is  simply  the  classical  trajectory  as  a  function  of 
the  initial  conditions  and  time. 

Moreover,  we  have  demonstrated  in  I  that 

a*(x,p\t)  =  oc(x,p;  t)  +  0(A2).  (14) 

Thus,  to  obtain  the  0(h)  corrections  to  an  expectation  value,  one  may  replace  a*  by  oc  in 
Eq.  (12).  Tb*«  a  substantial  simplification,  since  oc  cam  be  obtained  by  solving  Newton’s 
equations,  while  o*  would  require  the  solution  of  Heisenberg’s  equations  of  motion.  There¬ 
fore,  up  to  quantum  mechanics  enters  <  A(t)  >  only  through  the  initial  Wigner 

function  fq. 

In  I,  we  illustrated  the  use  of  Eq.  (14)  with  a  simple  one-dimensional  example.  In  the 
following  section,  we  apply  it  to  obtain  formulae  for  the  quantum  corrections,  A/j.,  A E^, 
and 

IV.  Quantum  Correction  Formulae 

For  incidon  energies  ~  1  ev,  the  Born-Oppenheimer  approximation  provides  an  accu- 
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rate  description  of  the  atom-surface  system.  We  thus  take  the  Hamiltonian  to  be 

2 

»  =  I); +  «'(*.{.')  + 5 1>? +  "?«?)•  (15) 

t 

where  p  =  (p*,Py,Pz)  is  the  incidon’s  momentum,  x  =  (z,y,z)  is  the  incidon’s  position, 
and  m  is  the  incidon’s  mass.  tt,  and  £,•  are  coordinates  for  the  normal  modes  of  the 
target  atoms,  whose  corresponding  quantum  operators  obey  the  commutation  relation 
[£«>*;]  =  * ui  is  the  frequency  of  the  *th  normal  mode,  and  V (x,  £,)  is  the  incidon- 
target  interaction  potential.11 

We  take  the  initial  state  of  the  incidon  to  be  a  Gaussian  wave  packet 

0(x)  =  0  *  (16) 


where  xq  is  the  initial  average  position  and  po  is  the  initial  average  momentum.  Eventually, 
we  shall  take  the  limit  -7  -*  0  so  that  rp  approaches  a  plane  wave.  From  Eq.  (13),  we  obtain 
the  Wigner  function  corresponding  to  rp, 


/(X'P) 


1  -[7(x-x0)3+;4-(p-po)2] 

(xh)3 


(17) 


The  initial  state  of  the  target  for  zero  temperature  is 


(18) 


and  its  Wigner  function  is 


1  *? 

/<(&.»() -n  3*  ■ 


(19) 


The  case  of  finite  temperature  is  discussed  in  Appendix  A. 


The  total  initial  wave  function  is  simply 


*(x,£)=«x)  **(<>). 


(20) 
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and  the  total  initial  Wigner  function  is 


^(x.p, &,*,-)  =/(x,p)  Mti'Xi). 


(21) 


We  now  consider  any  operator  A(t)  =  <J[p(t)],  depending  only  the  incidon’s  momentum 
p(t),  where  Q(p)  is  a  smooth  function  of  p.  We  apply  Eq.  (12)  to  Q[p(<)],  using  the  Wigner 
function  (21).  Exploiting  Eq.  (14),  we  expand  in  powers  of  tl  and  find 


<  <?[p(t)]  >=  Q)  *  J  dZX 

X  |gc(x,P0,O,O;t)  +  ^[^£  <  £7(0)  >  <*?( 0)  >]|  +0(h2), 


(22) 


where  Qc  is  the  classical  function  corresponding  to  Q[p(t)]  and  with  the  partial  derivatives 
being  evaluated  at  p  =  pg,  =  0,  and  irf-  =  0.  The  initial  mean  square  values  of  the 
normal  mode  coordinates  are  given  by 


<  £?(o)  >= 


2u /,’ 


<  *ii°)  >  =  2^* 

and  so  Eq.  (22)  gives  the  0(fl)  correction  to  <  Q[p(t)|  >.12  We  note  that 


(23) 


<?c(x,p,  £,-,*■,;*)  =  Q[Pc(x,p,£,-,ti;0),  (24) 

where  pc  is  the  classical  momentum  trajectory,  with  its  arguments,  (x,  p,  £,,  jt,-),  specifying 
the  initial  conditions. 


Defining 


=  (j)  ’  f  7(x  Xo)I<Jc(x,p0, 


(25) 


Eq.  (22)  can  be  written  in  the  compact  form 
<  9{p(*)]  >=Qc(x0,Po,0,0;t) 

l^\d2Qe 


+ 


a2Qc 


jEl-^r  <  «?(<>)  >  <  »?(o)  >]  +o{K 

i  'i  « 


(26) 
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We  assume  that  the  interaction  between  the  incidon  and  the  target  is  negligible  initially 
and  also  long  after  the  collision.  Thus,  as  t  — ►  oo,  Qc  becomes  independent  of  f,  since  the 
incidon’s  momentum  distribution  approaches  its  asymptotic  limit.  We  now  take  the  limit 
20  — *  oo,  ^  ►  0,  and  -72$  — »  00.  This  means  that  initially  the  wave  packet  is  very  distant 

from  the  target,  very  spread  out,  but  has  negligible  overlap  with  the  target.  In  this  limit, 
Eq.  (26)  is  evaluated  by  calculating  classical  trajectories  and  averaging  uniformly  over 
their  impact  parameters.  The  evaluation  of  partial  derivatives  such  as  those  appearing  in 
Eq.  (26)  is  discussed  in  I. 

We  now  use  Eq.  (26)  to  obtain  expressions  for  A /*,  A 2?*,  and  The  probability 

of  finding  the  incidon  with  a  momentum  p  at  a  time  t  is 

G(p)  =<  i3(p-p(t)]  >  .  (27) 


As  discussed  in  section  I,  G( p)  does  not  have  an  expansion  in  powers  of  h.  However,  J* 
and  A/n,  which  do  have  expansions,  can  be  expressed  in  terms  of  integrals  of  G. 

Define 

<Mp;M)  h  (£)"££  (pf  «3(p'  - 

(28) 

where 

p±  =  [2m(Ek±())K  (29) 

and  e  is  an  A-independent  energy  interval  <  |  —  Ey  |. 

The  smoothed  angular  distribution  of  Eq.  (5)  is  then 

7t(M)=<<?o[p(0;M>-  (30) 


and  the  energy  moments  of  Eq.  (8)  are  simply 


T7  -  <Qn(m;*,*)> 

<Qo(p(0;M> 


(31) 
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Applying  Eq.  (26)  to  Eq.  (30),  we  have 

?&,+)  =  7*,o(*,*)  +  i£[ 


-  T2  <  f?(°)  >  +  <  *?{»)  >]  +  0(A2),  (32) 


where  the  classical  smoothed  intensity  7*0  is  considered  as  a  function  of  the  initial  values 
of  the  normal  coordinates.  Taking  the  limit  a  -*•  0,  we  obtain  for  the  leading  correction  to 
the  angular  distribution 


^  =  i£[^<e7(o)>+^<*?(o)>]. 


(33) 


d(*  '  * r? 

Note  that  this  result  is  independent  of  the  energy  interval  c.  In  Appendix  B,  we  give  the 
explicit  expression  for  the  classical  intensity,  J*  0,  in  terms  of  classical  trajectories. 

Eq.  (33)  can  be  written  in  a  more  useful  form  by  introducing  coordinates  Xj,  represent¬ 
ing  physical  displacements  of  the  target  atoms  from  their  classical  equilibrium  positions, 
and  their  corresponding  momenta  Pj.  These  variables  are  linearly  related  to  the  normal 
coordinates  and  lead  to  the  expression 

A/*  -  5  E  >  +S&  <  *mm  >]  •  (34) 

i,j  J  J 

This  equation  formally  contains  an  infinite  number  of  terms.  However,  in  practice  often 
just  those  containing  derivatives  with  respect  to  positions  (or  the  corresponding  momenta) 
of  the  target  atoms  near  the  classical  impact  point(s)  contribute  significantly.  Hence,  only 
a  small  number  of  terms  is  needed,  making  the  evaluation  of  (34)  numerically  feasible. 

Expressions  similar  to  (34)  can  be  found  for  AEk  and  62Ek.  Applying  Eq.  (26)  to  Eq. 
(31)  we  obtain,  after  some  simple  algebra,  the  quantum  correction  to  the  position  of  an 
energy  peak 

I  d2Ek  d\n{Ik,o)  dEk] 


«  -gdisfsfc  *  - 


(35) 
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Similarly,  for  the  square  width  of  a  peak,  we  have 


c2El  y^rd-Et  dEk  y  (c\\Y  /ni  ^  ,  dEk  dEk 

k  ^  Lax,  ax,  Xx^x>  ^  >  +  a  Pi  dp , 


<  Pi(0)Pj( o) 


(36) 


As  with  Eq.  (34),  only  a  few  terms  of  the  sums  in  (35)  and  (36)  need  to  be  calculated  in 
practice. 


V.  Concluding  Remarks 

A  naive  view  of  a  crystal’s  quantum  mechanical  ground  state  is  that  it  is  similar  to 
a  classical  crystal,  but  that  the  normal  coordinates  fluctuate  randomly  with  a  Gaussian 
distribution  whose  widths  are  determined  by  the  second  moments  of  Eq.  (23).  In  fact, 
the  classical  phase  space  distribution  function  for  such  a  model  is  identical  to  the  Wigner 
function  of  Eq.  (19).  In  this  paper,  we  have  explicitly  shown  that  up  to  and  including  terms 
of  O(h)  the  quantum  corrections  to  the  incidon’s  final  angular  and  energy  distributions, 
A  It,  A  Ej,.,  and  are  correctly  given  by  this  picture. 

This  observation  allows  one  to  derive  simple  and  explicit  expressions  for  the  leading 
quantum  corrections  to  the  angular  distribution  [Eq.  (34)]  and  to  the  energy  distributi  n 
[Eqs.  (35)  and  (36)].  These  exact  expressions  depend  only  on  classical  trajectories  am' 
equilibrium  correlation  functions  of  the  target.  Our  results  have  practical  applications  to 
the  scattering  of  ~  1  ev  atoms,  for  which  a  semiclassical  approximation  is  appropriate.  A 
paper  describing  such  applications  to  the  scattering  of  rare  gas  atoms  from  alkali-halide 
surfaces  is  in  preparation. 
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Appendix  A:  Finite  Temperature 


In  this  appendix  we  extend  Eqs.  (34),  (35),  and  (36)  to  finite  temperature. 

If  the  target  is  at  a  temperature  T  >  0,  it  is  described  by  the  density  matrix  operator 

=  (A1) 


where 


A  =5  £(*? +*?(?), 


(A2) 


Z  =Tr(4,)> 


and  k  is  Boltzmann’s  constant. 


/  fiUi  \ 

r  ,  /  hu/v  \ 

“pr*“h(»r) 

(  h 

+  hw,)] 

The  Wigner  function  corresponding  to  pt>  defined  by  Eq.  (13)  with  the  projection 
operator  |t/>o  ><  r/>o|  replaced  by  pt,  is10 

.-v  m  j  •.  v  /(.I.  C2  v  S 

(A3) 

Fluctuations  in  the  positions  and  momenta  of  the  target  atoms  are  now  due  to  both 
quantum  and  thermal  effects.  The  leading  effect  of  these  fluctuations  on  the  expectation 
value  of  an  operator  (J(p(t)j  is  still  described  by  Eq.  (26),  provided  the  mean  square 
deviations  of  the  normal  coordinates  given  by  (23)  are  replaced  by  their  finite  temperature 
values 


<*?«,)  >=^co«h^. 


(A4) 


This  may  be  justified  by  repeating  the  analysis  of  section  IV.  Therefore,  Eqs.  (34),  (35), 

A  A  A  A 

and  (36)  are  still  valid  except  that  the  correlation  functions  <  X,X;  >  and  <  P, Pj  >  are 
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to  be  evaluated  for  the  temperature  T  instead  of  for  the  ground  state.  If  kT  is  not  much 
larger  than  hup,  they  correctly  give  the  leading  quantum  and  thermal  effects.  If,  however, 
kT  is  much  larger  than  hupi  then  Eqs.  (34),  (35),  and  (36)  give  the  leading  thermal 
effect,13  but  the  quantum  correction,  now  being  of  0(h2),  is  not  the  only  correction  of 
that  order. 


Appendix  B:  Classical  Intensity 


Here  we  give  an  explicit  expression  for  the  classical  intensity  Ik  q. 

Trajectories  with  a  particular  final  direction  correspond  (for  a  periodic  surface)  to  a 
finite  number  of  physically  distinct  impact  parameters  labeled  by  the  index  k.  Let  nk  be  the 
density  in  the  zy-plane  (number  per  unit  area)  of  a  particular  type  of  impact  parameters. 


For  a  given  initial  incidon  momentum,  the  classical  scattering  angles  are  functions 
0c{a,b)  and  <pc(a,b)  of  the  impact  parameters  (o,  b).  Defining 


D(a,b)  = 


m 


we  simply  have 

_  nk  D~l{ak,bk) 

*’°  s\n[Oc{ak,bk)]  ' 

where  ( aktbk )  are  the  impact  parameters  for  trajectory  k. 


m 
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FIGURE  CAPTIONS 


1)  The  xyz-coordinate  system.  The  surface  is  parallel  to  the  zy-plane.  The  point  A  on 
the  xy-plane  represents  an  impact  point  with  impact  parameters  (a,  6).  An  incidon  which 
begins  at  the  point  B  will  scattered  ^identically  to  an  incidon  of  the  same  momentum  that 
begins  at  B' . 

2)  Two  physically  distinct  trajectories  that  scatter  in  the  same  direction  but  with  dif¬ 
ferent  final  energies. 

3) The  intensity,  1(0, <j>;  E),  as  a  function  of  the  energy  E.  Classically,  it  consists  of  a 
sum  of  delta  function  peaks  (two  are  shown),  while  quantum  mechanically  these  peaks  are 
broadened  to  a  finite  width. 
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